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D , Abstract. Thermal photons can provide information on the temperature 

of the new state of matter created at RHIC. In the px region of 1— 
3 GeV/c thermal photons are expected to be the dominant direct photon 
source. Therefore, a possible excess compared to a pure decay photon sig- 
nal due to a thermal photon contribution should be seen in the double ratio 
(7/7(7r°))Mca S urcd/(7/7( 7r °))simuiated, if sufficient accuracy can be reached. We 
present a method to reconstruct direct photons by measuring e + e _ -pairs from 
iy~j ' external photon conversions. 
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1. Introduction 

Direct photons are produced during all stages of heavy ion collisions at the Rel- 
^ ■ ativistic Heavy Ion Collider (RHIC). Because they do not interact strongly, they 

escape the medium unaffected by final state interactions and provide a promising 
signature of the earliest and hottest stage of the quark-gluon plasma (QGP) [[I]. 

On a microscopic level, the main sources of direct photons from a QGP are 
quark-gluon Compton scattering (qg — > 75), quark- antiquark annihilation (qq — > 
717) and bremsstrahlung involving thcrmalized partons [ 2 . Direct photons are also 
produced in initial hard scattering processes which involve the same reactions but 
among the incoming partons. 

At RHIC energies thermal photons are predicted to be the dominant source of 
direct photons in a px window between 1-3 GeV/c [Q]. 

Direct photons have been measured with PHENIX in Au + Au collisions at 
t/s~n~n = 200 GeV [ [3] . The inclusive photon spectra measured with the Electro- 
magnetic Calorimeter (EMC) have been compared to the expected background from 
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Fig. 1. Double ratio of the measured invariant yield ratio, (7/^°). . ,, to the 

° J ' V " / Measured' 

background decay ratio, (7/ 7r °) BackgI . ound ; as a function of pr for minimum bias 
and for five centralities of Au + Au collisions at ^/s^n = 200 GeV. 



hadronic sources, based on the measured ir° and r\ spectra and a cocktail of other 
hadronic decays (77', Kg, uj), assuming ttit scaling. 

Fig. [T] shows the double ratio of the measured invariant yield ratio to the back- 
ground decay ratio as a function of pr for minimum bias and for five centrality 
classes. The measurement of direct photons production at high px scales with the 
number of binary collisions in agreement with NLO perturbative QCD predictions 
and therefore confirms medium effects as the origin of jet suppression. In the low pt 
region, where a thermal signature is expected, a significant measurement remains 
limited by systematic uncertainties due to the energy resolution and the photon 
identification with the EMC. 

In order to overcome such limitations, dielectron pairs offer some advantages 
because of the superior momentum resolution of charged particles at low momenta 
and excellent identification of conversion photons; while other methods [3] try to 
use low mass dielectron pairs from internal conversions, the method presented here 
uses real photon conversions in the beam pipe. 
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2. Thermal photon analysis 

The excellent capabilities of the PHENIX detector to measure electrons suggest to 
circumvent the limitations of the conventional direct photon measurement [[3j at low 
photon energies by measuring photons via their conversion into e + e _ -pairs. The 
momentum resolution (~ 1 %) of the charged tracking devices proves superior to 
the energy resolution of the EMC (~ 10 %) in the pr region of interest (1-3 GeV/c). 

To identify e + e _ -pairs from photon conversions, a single electron identification 
cut is applied, which require signals from at least two phototubes in the Ring 
Imaging Cherenkov Detector (RICH) matching to a reconstructed charged track in 
the Drift Chamber (DC). No further electron identification cuts were applied since 
the pair cuts (see Sect. I2.1|) to separate conversion photons from other e + e~-pairs 
are more efficient and powerful enough to provide a very clean photon conversion 
sample. 

The extracted photon conversions are tagged with photons reconstructed in the 
EMC to determine the contribution from ir° — > 77 decays (see Sect. I2.2|) . 

All yields are measured as a function of pt of the e + e~-pair, which makes a 
direct comparison of the inclusive photon yield, A/l ncl , and the tagged photon yield, 

A^° tag , possible: 

e e+e - a e+e - 7 mcl (PT) (1) 
e e+e - a e+e - e 7 (jp T ) f 7^° (Pr) (2) 

The measured yield of inclusive photons depends on the reconstruction efficiency 
e e + e - and the PHENIX acceptance a e + e - of the conversion e + e~-pair . The tagged 
photon yield depends in addition on the efficiency to reconstruct the second photon 
in the EMC e 7 (pr) and on the conditional probability / to find it in the EMC 
acceptance, given that the e + e~-pair has been reconstructed already. Here, e 7 (py) 
is weighted with the px distribution of the e + e~-pair. In the ratio 7V 7 ™ c/ /A^J r °* a9 
the e + e~-pair reconstruction efficiency and acceptance correction factor cancel. 

A ratio of the hadronic decay photon yield, N^ dT , and the tagged photon yield 
from 7T° decays, 2V7 tag , is calculated with simulations. 

Nf^( PT ) = fNf(p T ) (3) 

The comparison of the ratio in data and in simulations in a double ratio leads 
to an expression that is equivalent to the ratio of inclusive and decay photons as 
shown in Eq. ((4]). 

7 hadr (PT) / N^( PT ) \ 

I 7 Nf{p T ) J 

\ ' / Sim 



Nf tas ( PT ) = 
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Fig. 2. Invariant mass of e + e -pairs before (black) and after (red) applying cuts 
on the orientation of the e + e _ -pair in the magnetic field. 



The only remaining factors are the reconstruction efficiency of the photon in the 
EMC, e 7 (pr)j and the conditional acceptance / in the simulation part of the dou- 
ble ratio, which have both been determined with Monte Carlo simulations (see 
Sect. EH). 

2.1. Photon Conversions 

Since the PHENIX tracking algorithm assumes the track to originate from the col- 
lision vertex, off-vertex conversion pairs are reconstructed with an artificial opening 
angle which leads to an invariant mass that is proportional to the radius at which 
the conversion occurs. 

Therefore, photon conversions that occur in the beam pipe material (Be, 0.3 % 
radiation length) at a radius of 4 cm are reconstructed with an invariant mass of 
~ 20 MeV/c 2 . Fig. [2] shows an invariant mass spectrum of e + e~-pairs in the range 
0-0.1 GeV/c 2 . The peak from photon conversions in the beam pipe at 20 MeV/c 2 
can be clearly separated from Dalitz decays ir° — ► 7e + e~, which dominate the spec- 
trum below 10 MeV/c 2 , and combinatorial background pairs, whose contribution 
increases toward higher invariant masses. 

The photon conversion pairs, which have no intrinsic opening angle, can be 
distinguished from Dalitz decays and purely combinatorial pairs by cutting on the 
orientation of the e + e _ -pair in the magnetic field. 

Fig. [5] shows the invariant mass spectra of e + e~-pairs before (black) and af- 
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ter (red) applying these pair cuts. The yield from integrating the mass region 
< 35 MeV/c of the conversion peak is corrected for the remaining px dependent 
contamination of ~ 15.0 ± 2.0 (syst) % due to combinatorial e + e~-pairs which has 
been determined with mixed events. 

2.2. Tagging of Decay Photons 

To reveal which of these conversion photons come from 7r° — ► 77 decays, the e + e _ - 
pairs in the conversion peak are combined with photons which have been measured 
in the EMC, under loose cuts based on the time of flight and the shower profile for 
photons with a minimum px of 0.3 GeV/c, and their invariant mass is calculated 
(see Fig.©. 

The reconstruction efficiency e 7 (pr) of the loose photon has been estimated 
with a full GEANT simulation which embed simulated photons into real EMC 
data, therefore providing a combined information on the photon identification effi- 
ciency and occupancy effects. The overall efficiency is determined to be 82 ± 1 % 
independent of px beyond the minimum px cut off. 

Conversion photons that are identified as decay products of ir° can be tagged as 
iV* tag . This signal has a large combinatorial background due to the high photon 
multiplicity in Au + Au collisions. 

The combinatorial background is reproduced with an event mixing method, 
which creates uncorrelated pairs of photons and e + e _ -pairs from different events. 
The mixed event spectrum is normalized to the same event spectrum well outside 
the 7T° mass region (0-100 MeV/c 2 , 170-250 MeV/c 2 ) and subtracted. 

The statistical error on the normalization factor is on the order of 0.2 % and 
depends only on the statistics in the same event spectrum in the normalization 
region. As an example, the resulting ir° signal for e + e _ -pairs with 0.8 < px < 
1.2 GeV/c is shown as insert in Fig. [3] 

Mean and er are determined by a fit of the background subtracted data with 
a Gaussian. The data are also fitted to the sum of a second order polynomial and 
a Gaussian, to take into account the possibility that the shape is not completely 
described by the mixed event spectrum. The difference in the resulting mean and a 
is negligible. The mean and a obtained by the fit are then used to integrate the data 
in a region ±1.5 a around the mean, chosen to optimize the signal to background 
ratio. 

The statistical error on the extracted tt° signal is given by: 



With FG(i) and BG'(i) being the yields in bin i of invariant mass spectrum in same 
events and normalized mixed events, respectively, the summations are performed 
over the integration region. It is important to note that the last term in Eq. (|5j). 
is the square of the sum over the normalized background, and therefore, depends 
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ye + e" inv. mass with: 0.8 GeV < e + e" pair-p T < 5.0 GeV 




x'/ndf 25.32/22 
Yield 9300± 475.7 

Mean 0.1 382 ± 0.0009 
Sigma 0.01 51 6 ± 0.00087 



0,25 
m ye . e . [GeV/c] 



Fig. 3. Invariant mass of 7e + e _ -triplets in same events (black) and normalized 
mixed events (red) for e + e~-pairs with 0.8 < pr < 5.0 GeV/c. The insert shows 
the invariant mass of 7e + e~~triplets after background subtraction for e + e~-pairs 
with 0.8 < pt < 1-2 GeV/c. A fit with a Gaussian is drawn and the resulting 
parameters shown in the box in the upper right of the graph. 



on the integration region and is not bin independent. Different integration regions 
have been used. Variations in the resulting yield have been used to set a systematic 
uncertainty on the yield extraction of 2.5% independent of px- 

The loss of tag due to the external conversion of the second photon is cor- 
rected by a factor 1 — p con v = 94 ± 2 %. In this factor p con v is the conversion 
probability due to the material budget between the vertex and the Pad Chamber 3 
(PC3) in front of the EMC. 

2.3. Simulations 

The contribution of hadronic decays has been determined with a fast Monte Carlo 
simulation of ir° and r\ Dalitz decays. A parameterization of the ir° spectrum 
measured by PHENIX [ [5] has been used as input. The r\ distribution has been 

generated assuming mj scaling [px — > \jp\ + m r) ~ m Jo) of the tt° spectral shape 

and a normalization at high px to rj/ir = 0.45 ± 0.04, according to PHENIX data [ 
IH1 E]. The relative error of 9 % on the rj/ir° ratio is reduced by the branching ratio 
of the two photon decay and results in a 3 % error in the ratio A^ adr /iV^ r . 

The contamination due to neutral Kaons which decay before the beam pipe has 



Measurement of photons via conversion pairs with PHENIX at RHIC 



7 



J 2.4 
£ 2.2 
2 
1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 



— Direct Photons in MinBias Au+Au collisions at\ls = 200GeV 


L PHENIX Preliminary 


• Tagging method 




• Statistical method 




^External Conversions 


| I 

-, i , , , , i , , , , i , , , , i , , , , 


• 

t • 


, , , , i , , , , i , , , , i , , , , i 



1.5 



2.5 



3.5 



4 5 5 
p T [GeV/c] 



Fig. 4. Ratio of inclusive photons to decay photons, 7 incl {pt) /7 hadr (pt), in Min- 
imum Bias Au + Au collisions at ^/snn — 200 GeV, for three different methods 
methods. 



been found negligible (~ 1 %) and has been folded into the systematic error on the 
simulations. 

The conditional probability / that the photon is reconstructed in the EMC 
once the e + e~-pair is reconstructed already was calculated with a fast Monte Carlo 
simulation of 7r° — > je + e~. The use of Dalitz decays is justified by the fact that 
the pr spectra of photons from 7r° — > 77 are essentially identical to the e + e~-pair 
Pt spectrum from 7r° — > je + e~ . 

After the e + e~-pair has been filtered in the detector acceptance, the conditional 
acceptance / for the second photon has been calculated taking dead areas of the 
detector into account. Uncertainties in calculating / are found to be 5 %, which is 
the largest source of systematic errors. The limited energy resolution of the EMC 
of <?e/E = 5 % © 9 %/ \f~E introduced an additional systematic error of ~ 1 % due 
to the pt cut at 0.3 GeV/c. 

3. Conclusions 

Fig. |4] shows a preliminary result for the double ratio 7 incl (pr) /j hadT {pt) as in 
Eq. |4]for Minimum Bias Au + Au collisions at ^/snn — 200 GeV. The main sources 
of systematic errors arise from the uncertainties in the description of the detector 
active areas, in the peak extraction and in the assumptions of the ir° shape and 
give a final systematic error on the double ratio of ~ 7%. 
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The result of the presented analysis is compared to two other direct photon 
measurements in PHENIX. The first one [ [5] is based on the same tagging method, 
but instead of photons coming from conversions in the beam pipe, the clean photon 
sample is determined by selecting EMC clusters with very strict photon identifica- 
tion cuts. The second one is the ratio obtained from the statistical subtraction of 
measured ir° pt spectra from the measured inclusive photons [[9]. 

All three analyses indicate an excess due to direct photons above the decay 
photon spectrum in the px region 1-3 GeV/c. Despite of the statistical limitations 
due to the low conversion probability of 0.2 % this analysis offers a smaller system- 
atic error with respect to the conventional methods, and therefore, the intriguing 
ability to extend the extraction of a significant direct photon yield to low p?. 
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